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1-justification du choix de projet : 


Le ciment est un matériau essentiel pour la construction, qui est un des 
premiers secteurs d'activité au monde et un des employeurs les plus 
importants. D'’importance vitale pour le logement et les infrastructures de 
base, la filière ciment joue un rôle clé dans le développement économique et la 
réduction de la pauvreté des pays émergents. Cependant, l'industrie 
cimentière est aussi un des plus gros émetteurs de dioxyde de carbone et a 
d'importants impacts sociaux et environnementaux. Très capitalistique, elle 
requiert de lourds investissements qui nécessitent des financements et une 
rentabilité sur le long terme. Enfin, ce secteur, très énergivore, dépend des 
cycles de l'économie et de la santé du secteur de la construction ; cela rend les 
charges d'exploitation et le chiffre d'affaires très volatiles. Dans ce contexte, 
les institutions financières de développement (IFD) jouent un rôle crucial dans 
le soutien aux projets cimentiers sur les marchés émergents. Par leur 
participation et leurs conditions d'accès au financement, elles peuvent rendre 
les projets cimentiers plus durables et favoriser les innovations. 

L'industrie du ciment est très capitalistique et exige des financements à long 
terme auxquels il est difficile d'accéder dans les pays en développement. 

En favorisant donc la disponibilité des matériaux de construction, les IFD 
participent au dynamisme du secteur. La capacité des pays en développement 
de construire les infrastructures et les logements dont ils ont besoin pour 
réduire la pauvreté et soutenir la croissance dépend, à son tour, de la bonne 
santé du secteur de la construction. Enfin, de meilleures infrastructures 
conduisent à la croissance du PIB, créent des emplois, renforcent les PME et 
encouragent les autres investisseurs étrangers. 


D'où l'intérêt du développement du processus de production 
cimentière et son optimisation. 
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2- Introduction à l’industrie cimentièere : 


Le béton, le matériau de construction le plus utilisé au monde, est composé de 
60 à 80% de granulats et de 20 à 40% de liant (ciment Portland et eau). C'est 
Louis-Joseph Vicat qui a été le premier à fabriquer artificiellement de la chaux, 
mais c'est Joseph Aspdin qui a breveté le ciment Portland en faisant chauer un 
mélange d'argile et de calcaire. Au début des années 2000, on comptait, au 
Canada, 16 cimenteries qui fabriquaient ensemble 12,8 million de tonnes de 
ciment annuellement. Le ciment est composé essentiellement d'oxydes de 
calcium, de fer, d'aluminium (alumine), de silicium (silice) et de soufre 
(sulfate). Ces oxydes proviennent de matières premières comme le calcaire 
(source de chaux), l'argile (source de fer, de silice et d'alumine) et le gypse 
(source de sulfate). Certains de ces matériaux, dont le calcaire et l'argile, 
peuvent être extraits d'une carrière à proximité de la cimenterie. 


Calcaire + Argile 
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Deux types de procédé existent pour fabriquer le ciment Portland, le procédé 
par voie humide et le procédé par voie sèche. Aujourd'hui, la grande majorité 
des cimenteries utilisent le procédé par voie sèche. Dans ce procédé, la pierre 
provenant de la carrière est d'abord concassée puis envoyée dans un silo de 
pré-homogénéisation. A la sortie de ce silo, cette pierre est dosée avec les 
autres matières premières an de produire un mélange avec la composition 
désirée. Ces matières premières sont ensuite pulvérisées et mélangées. Le 
mélange ainsi obtenu est appelé mélange du cru de ciment, est composé 
essentiellement d'oxydes de calcium, de silicium, de fer et d'aluminium et est 
homogénéisé dans des silos. Le cru est envoyé à un four rotatif où se produit 
une calcination qui transforme chimiquement le cru en clinker de ciment. Il 
reste ensuite à pulvériser le clinker avec du gypse pour obtenir un ciment 
Portland prêt à l'utilisation. 
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3- Le but du contrôle du cru de ciment : 


Le but de ce contrôle est d'obtenir une entrée au four qui produira un clinker 
de ciment de qualité tout en minimisant la consommation d'énergie du four. 
Etant donné que la composition du clinker doit atteindre des cibles strictes, 
c'est aussi le cas du mélange de cru. Le dosage des différents matériaux à 
entrer dans le mélange est contrôlé par la vitesse des courroies qui 
transportent ces matériaux. Un instrument d'analyse de la composition du 
mélange est placé soit à l'entrée ou à la sortie du broyeur et permet d'obtenir 
cette composition sur une base régulière. Les nouveaux instruments d'analyse 
du mélange tels que les instruments aux rayons gamma (prompt gamma 
neutron activation analysis ou PGNAA) permettent une fréquence 
d'échantillonnage plus élevée (de l'ordre de quelques minutes plutôt que 
quelques heures). 


4- Notions de base de l'industrie du ciment : 
Le clinker est composé essentiellement de quatre minéraux. La composition du 
clinker, et donc du cru de ciment, est souvent exprimée en fonction de ces 


quatre phases, 





Minéral Abréviation 
































silicate tricalcique (alite) C3S 
silicate dicalcique (bélite) C2S 

aluminate dicalcique C3A 
aluminoferrite tétracalcique C4AF 








Les équations de Bogue permettent de calculer le pourcentage en poids de 
chacun des minéraux qui sera présent dans le clinker à partir de la composition 
en oxydes de calcium (CaO ou C), de silicium (SiO2 ou S), d'aluminium (AÏ203 
ou A) et de fer (Fe203 ou F) et sont résumées dans la matrice de l'équation. 
Ces équations ont été déterminées en présumant que tout l'oxyde de fer se 
combinera avec les autres oxydes pour former l'aluminoferrite té- 15 
tracalcique, que toute l'alumine restante se combinera avec la chaux pour 
former l'aluminate dicalcique, que toute la silice se combinera avec la chaux 
pour former la bélite et qu'enfin le surplus de chaux se combinera avec la 
bélite pour former l'alite. 


Cas | 4.017 T.60 1, 43 
C,5 


C3A 0,00 0,00 —1.69 


JE E LO 


3,07 8,60 1,0% 5, 07 


— 


en 
ER 
b 7 i 
TERME PR ES RC RE | 


CaAF 0, 00 0, CN) J LM 0, CM) 


Il est nécessaire de définir les vecteurs suivants : le vecteur contenant les 
compositions du mélange, xm et les vecteurs contenant les compositions des 
différentes courroies, xi , où i représente le numéro de la courroie. pour i = 
1...nc et nc est le nombre de courroies. 


Dans cette étude, l'analyse du cru se fait avant le broyeur. D'ailleurs, Berg 
(1999) a noté qu'en faisant l'analyse avant le broyeur, les délais entre le point 
de contrôle et le point d'analyse étaient plus faibles que si l'analyse était faite 
après le broyeur, ce qui permet d'améliorer le contrôle. Weiss et Gartner 
(1986) ont aussi noté une amélioration du contrôle lorsque l'analyse était 
effectuée avant le broyeur sur des simulations. Afinn de simplifier l'étude, ces 
délais ont été complètement négligés. Le mélange simplifié est donc modélisé 
comme suit : 


xm = XF où F=] f1...fi...fnc]| avec fi représente la proportion de la 
ième courroie et X =[ x1...xi...xnc] 

xi est la composition de la ième courroie telle que présentée à l'équation.Afin 
d'assurer que la somme des proportions des courroies donne 100% 

F'1 = 1 avec le vecteur 1 de taille nc x 1. 


5- Optimisation : 


L'objectif du contrôle du mélange de ciment au cru étant d'avoir un matériau le 
plus homogène possible, le contrôle se fera en minimisant l'erreur au carré par 
rapport aux cibles posées pour les composantes du mélange. L'optimisation 
doit nous permettre d'obtenir les proportions des courroies donnant un 
mélange dont les composantes principales se rapprochent le plus des cibles. 
Etant donné que les proportions de quatre courroies sont recherchées par 
l'optimisation, nous pouvons chercher à minimiser l'erreur au carré par rapport 
à 3 ou 4 cibles tout en contraignant la somme des proportions des courroies à 
être égale à 1. Si l'optimisation se fait sur 4 cibles, la fonction d'optimisation 
s'écrit : 

J(xm(t}) = (xm(t) — x0)(xm(t) — xO)t 

(xm(t) - x0)t : la transposée de (xm(t) - x0) 


Si l'optimisation se fait sur 3 cibles, la fonction d'optimisation est semblable, 
mais xm et x0 ne comportent plus que trois des quatre oxydes initiaux. De 
plus, les oxydes initiaux peuvent être remplacés par 3 où 4 des constituants du 
ciment définis. Les valeurs du mélange optimal, x*m(t), sont fonction des 
proportions optimales, F* des différents matériaux arrivant au mélange. 
Lorsque les compositions des courroies au temps t sont connues, la relation 
entre 18 le mélange optimal et les proportions optimales est donnée par 
l'équation : x*m(t) = Xest(t) F*(t) où Xest(t) est la matrice des compositions 
des courroies qui peut être estimée par un filtre. 
Par contre, lorsque les compositions des matériaux au temps t ne sont pas 
connues, la relation entre le mélange optimal et les proportions optimales est 
donnée par l'équation. 

x$ (t) = x (E) — XË + XF" (t) 
où x estm (t) représente la composition estimée du mélange au temps t, 
X représente la matrice des compositions des courroies moyennes et F 
représente la moyenne des proportions des courroies utilisées dans 
l'estimation du mélange au prochain pas de temps. Les compositions des 
courroies et la composition du mélange au temps t peuvent être estimées à 
partir des données du mélange connues (du temps 1 au temps t-1) grâce à un 
filtre (estimateur) développé dans le domaine de la géostatistique et nommé 
cokrigeage. Lorsque l'équation est utilisée pour dans le calcul des proportions 


optimales, un gain, G, peut être appliqué sur l'optimisation. L'équation devient 
alors : 

xt (t) = xt (t) — GXF + XEF (t)) 
où G est une constante à déterminer . Ce gain permet d'appliquer seulement 
une portion de la correction. Enfin, la somme des proportions des courroies 
doit être égale à 1. 
Donc, pour nc courroies 


Xss 


i=] 


Pour des questions de simplicité, ceci se fait en posant que la proportion de la 
dernière courroie est égale à la différence entre 1 et la somme des proportions 
des autres courroies. 


Comme la méthode du simplexe de Nelder-Mead est utilisée pour 
l'optimisation, l'optimisation se fait sans contrainte. Par contre il est possible 
de rajouter des pénalités à la fonction d'optimisation qui sont liées à la 
variation des proportions des courroies et pour empêcher d'avoir des 
proportions des courroies négatives. La fonction d'optimisation devient donc 


J(Xmlt)) = (Xm(t) — x0) (Xml) — x0)" + où fneg + af var 


fneg = $  L(t)fi(E) où h(t) =1 si f(t) < 0 et 0 sinon 


i=] 


fvar = Y (F#4) — filt- D? 
i=l 


pour empêcher toute valeur négative, tandis que a2 est une valeur à 
déterminer. 
où ai et a2 sont des facteurs de pondération et habituellement ai = 1 x 109 


pour empêcher toute valeur négative, tandis que a2 est une valeur à 
déterminer. 


Dans les usines de ciment, un seuil minimum est souvent établi pour les 
proportions des courroies. Ainsi, si les proportions optimales sont plus faibles 
que les minima, la proportion de la courroie pourra soit être fixée au minimum 
soit être arrêtée (fixée à zéro). Dès que l'optimisation retourne la proportion 
d'une courroie plus faible que le seuil établi pour cette courroie, plusieurs 
scénarios doivent être pris en compte. Le premier scénario à l'étude consiste à 
fixer la valeur de la courroie qui était plus faible que le seuil à ce seuil. Ensuite 


le programme mathématique est minimisé sur les trois autres courroies. Le 
deuxième scénario consiste à fixer la valeur de la courroie qui était plus faible 
que le seuil à zéro et à minimiser le programme mathématique pour les trois 
autres courroies. 

J(xm(t)) = (xm(t) — x0)(xm(t) — x0) T + alfneg + a2fvar 
où ai et a2 sont des facteurs de pondération comme définis 


fneg — ; Tilti filt) Frère) F où Jilt} = 1 si filé) << O0 et O sinon 
=] 


fuar = 5 (fité] Filt Yy 
i— ] 


où Fmin, de taille nc x 1, représente les seuils sur chaque courroie. Les deux 
valeurs ainsi obtenues sont ensuite comparées an de choisir le scénario qui 
permet d'obtenir la valeur de J(xm(t)) minimale. Cette méthode qui consiste à 
travailler sur une seule courroie est une heuristique puisque de meilleurs 
résultats pourraient être obtenus en fixant une, ou plusieurs, autres courroies 
à leur seuil minimal ou à zéro. Elle a été choisie pour des raisons de rapidité de 
calcul. En effet, avant d'obtenir le minimum réel, il faudrait comparer les 
valeurs de : 


i possibilités lorsque l'on a nc courroies (et ( nc i ) représente le nombre de 
combinaisons possibles de i dans nc). Cette somme représente toutes les 
possibilités qu'une jusqu'à nc — 1 courroies soient fixées au seuil minimum ou 
à zéro, puis que la minimisation soit effectuée sur les autres courroies. 


Afin d'optimiser le programme mathématique le problème suivant doit être 
résolu. 
min F(t) (J(xm(t))) 


La méthode utilisée pour résoudre ce problème est le simplexe de Nelder- 
Mead. Cette méthode d'optimisation est une méthode de recherche directe. 
Elle est robuste puisqu'elle est tolérante au bruit dans la fonction à optimiser 
(Dennis et Woods, 1987). De plus, comme toute méthode de recherche 
directe, elle ne nécessite pas de connaître les dérivées. 

Dans l'algorithme de contrôle du ciment présenté ici, les conditions d'arrêt 
suivantes sont utilisées (lorsqu'une de ces conditions est respectée, 
l'algorithme d'optimisation s'arrête sur le simplexe de la dernière fonction 
évaluée) : tolérance sur la fonction à minimiser de 10—4 tolérance sur les 
variables de 10-4 nombre maximal d'itérations de 200 fois le nombre de 


variables. Et comme ici, la fonction a soit 3 ou 4 variables, le nombre maximal 
d'itérations est soit de 600 ou de 800. nombre maximal d'évaluations de la 
fonction de 200 fois le nombre de variables. Même si la méthode du simplexe 
de Nelder-Mead est une des méthodes de recherche directe les plus utilisées, il 
n'existe pas de résultat quant à la convergence globale de la méthode (Lewis 
et al., 2000). Par contre, cette méthode est disponible sur Matlab et permet 
d'incorporer facilement des contraintes de positivité des proportions dans la 
fonction objectif sous forme de pénalités 


6-ALGORITHME D'OPTIMISATION DE NELDER-MEAD : 


Il s'agit de générer n+1 sommets x1, x2,...,xn+1 du 
simplexe (dont nous pouvons calculer les valeurs f(x(1)), f(x(2)),..., f(x(n + 
1))). Ces points sont générés en utilisant l'estimation initiale du point minimal, 
x0 et en ajoutant 5% à chacune de ces composantes. Le simplexe initial est 
donc formé du point x0 et de ces n nouveaux vecteurs. 


:Il s'agit d'ordonner les sommets 
tel que f(x(1)) < F(x(2)) < ... < f(x(n + 1)), puis de déterminer le pire point 
x(n + 1). À chaque étape de l'itération, le pire point x(n + 1) est retiré au 
bénéfice d'un nouveau point accepté dans le simplexe. 


.Il s'agit de faire la réflexion du pire point, x(n + 1) déterminé à 
l'étape précédente par rapport au centroïde des autres points x(1}), x(2),..., 
x(n). Le centroïde, m, est donné par : ee 


La réflexion, r est ensuite calculée. r = (1 + p)m — px(n + 1) et p = 1 ( Puis, 
il faut calculer f(r). 


: Si f(x(1)) < f(r) < f(xn), alors le nouveau simplexe est 
déni comme étant < x(1)...,r,...,x(n) > et l'itération est terminée. 


:Si f(r) < f(x(1)), le sommet r peut être agrandi pour 
déterminer le point s. s = m + 2(m — x(n + 1)) et calculer f(s). 
(a) Si f(s) < f(r), alors s est accepté dans le nouveau simplexe. 
(b) Sinon, r est accepté dans le nouveau simplexe. 


:Si f(r) > f(x(n)), alors une contraction est nécessaire entre 
m et le meilleur des points entre x(n + 1) et r an d'obtenir le point c. 


(a) Sif(r) < f(x(n + 1)), alors la contraction est calculée à partir de 

c= m + (r — m)/2 et il faut calculer f(c). Si f(c) < f(r), le simplexe devient 
< X(1)...,c,...,x(n) > et terminer l'itération. 

Sinon, continuer à l'étape de rétrécissement. (b) si f(r) 2 f(x(n + 1)), la 
contraction est calculée à partir de cc = m + (x(n + 1) - m)/2 

Calculer f(cc) Si f(cc) < f(x(n+1)), alors le simplexe est 

< X(1).,cCc,...,x(n) >. Sinon, continuer à l'étape de rétrécissement. 


:Le rétrécissement du simplexe est calculé comme suit 
v(i) = x(1) + (x(i)— x(1))/2Vi=2...n+1 Et le nouveau simplexe devient 
< X(1)},V(2),...,vV(n+1)> 


: Une nouvelle itération peut commencer en retournant à l'étape de 
réflexion ou si un critère d'arrêt est respecté, l'algorithme est terminé. 


7-Lafarge Quality Mix Control software - logiciel de contrôle 
du mélange du cru du ciment de Lafarge au Maroc: 


Le Lafarge Quality Mix Control software (QMC) est le logiciel de contrôle du 
mélange du cru de ciment utilisé dans plusieurs cimenteries de Lafarge et 
auquel l'algorithme de contrôle du cru présenté dans ce mémoire sera 
comparé. Le fonctionnement de ce logiciel est présenté dans un document 
d'aide qmc (2004-2008). Ce logiciel sert à optimiser une fonction semblable à 
celle présentée en ce qui précède. 


f(xm(t)) = K(xm(t) — x0)(xm(t) - x0) T 


xm(t) est un vecteur 3 x 1 de 3 des oxydes, modules ou phases du mélange 
et x0 est aussi un vecteur 3 x 1 des valeurs cibles de ces oxydes, modules ou 
phases et K est une matrice diagonale contenant un coefficient, ki pour chaque 
matériau. 


Dans le logiciel QMC, une pénalité est aussi ajoutée à la fonction à minimiser, 
ce qui donne : 


f (mt) = K(km(t) — Xo)(Xmlt) — x0)” + kC FoP! 


C est un vecteur 4x1 contenant les coûts des quatre matériaux dans le 
vecteur F opt et k4 est un coefficient associé à la pénalité apportée par le coût 
des matériaux. Lorsque K est une matrice identité 3x3 et que k4 est nul. Dans 
ce cas, la minimisation utilisée pour obtenir les proportions des courroies 
optimales est définie comme en 3.16 Par contre, l'optimisation dans QMC se 
fait sous contraintes, tandis que l'optimisation dans l'algorithme développé 
pour ce mémoire se fait sans contraintes. En effet, l'algorithme d'optimisation 
de Nelder-Mead a été développé pour de l'optimisation sans contraintes. Les 


contraintes dénies pour l'optimisation dans QMC sont des contraintes 
d'inégalité sur les éléments chimiques mineurs du mélange : 


NasOogmin < N 20eg < Nan: maz 
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des contraintes d'inégalité sur les proportions des courroies 
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où n représente le nombre de 
courroies et dans le cas qui nous intéresse, n = 4. Dans QMC, la 
minimisation se fait aussi sous une contrainte d'égalité : 


fi 


* fi = 100 


i=] 


8-Comparaison en QMC utilisé en industrie cimentière au 
Maroc et l'algorithme de Nelder mead : 


L'optimisation dans QMC se fait soit grâce aux algorithmes d'optimisation de 
Harwell ou IMSL ou à un algorithme qui a été développé chez Lafarge nommé 
algorithme MA. L'opérateur peut choisir d'activer un ou plusieurs des 
algorithmes d'optimisation. Puis, le logiciel choisira l'optimisation en fonction 
de l'algorithme qui a retourné le meilleur résultat. La principale différence 
entre QMC et l'algorithme programmé dans le cadre de ce projet est la façon 
de calculer la composition du mélange au prochain pas de temps et la 
composition des courroies. 

Enfin, un gain existe dans le logiciel de contrôle QMC an d'atténuer les 
corrections apportées par le programme. Ce gain est à la discrétion de 
l'utilisateur. Dans QMC, le gain se définit comme suit : 

fi = fi(t-1) + G(f*i — fi(t-1)) où G représente le gain et f'i est la nouvelle 
proportion envoyée à la courroie. En d'autres mots, 

nouvelle proportion = ancienne proportion + gain (proportion optimale - 
ancienne proportion). 

Contrairement au gain présenté à l'équation ci avant, le gain utilisé dans QMC 
est à l'extérieur de l'optimisation. Ainsi, plus le gain est élevé, plus la 
correction sera élevée et plus le risque de corriger plus que nécessaire 
augmente tandis que plus le gain est faible, moins le risque de sur-corriger 
sera grand, mais la correction aura peu d'effet. Afin de comparer l'algorithme 


de contrôle développé dans le cadre de ce mémoire à QMC, un algorithme 
comparable à QMC a été programmé. Cet algorithme n'est que l'algorithme de 
contrôle présenté avant , mais dont quelques paramètres ont été choisis an de 
ressembler le plus possible à QMC. Donc l'optimisation se fait sur 3 cibles et le 
mélange est estimé à partir de l'équation où x est m est estimé à partir d'un 
cokrigeage ordinaire dont les données se résument à la dernière donnée de 
composition du mélange obtenue. Le gain présenté à l'équation 


Xmlt) = Xm (t) G(ŠË + XF (t) ) 


IYI 


est de 1 puisque ce type de gain n'est pas présent dans QMC. Enfin, étant 
donné que l'optimisation dans QMC se fait sous contraintes et que l'algorithme 
d'optimisation de Nelder-Mead ne permet pas de contraintes, les contraintes 
sont remplacées par des pénalités telles que présentées avant. La pénalité sur 
la variation des proportions des courroies a un coefficient, a2 nul puisque ce 
genre de pénalité n'est pas nécessairement présent dans QMC. Par contre, la 
pénalité sur les valeurs négatives des proportions des courroies a un coefficient 
de, a1 d'environ 109 an de respecter les contraintes sur les proportions des 
courroies présentes dans QMC. L'algorithme de contrôle avec les paramètres 
tels que dénis ci-haut équivaut approximativement au logiciel QMC lorsque, la 
matrice des coefficients, K de l'équation 3.50 est une matrice identité ; que le 
coefficient k4 de l'équation 3.52 est nul ; que les valeurs minimales pour les 
oxydes mineurs et les proportions des courroies dans les équations 3.53 et 
3.54 sont fixées à 0 tandis que les valeurs maximales sont fixées à 100 et que 
le gain, G, de l'équation 3.56 est fixé à 1. 


9-ÉTUDE D'UN CAS SYNTHÉTIQUE : 


:procédé de mélange de cru de ciment. 











Optimiseur p 





Ce procédé consiste en un mélange composé à partir des proportions des 
différents matériaux. Ces proportions ont été optimisées an d'obtenir un 
mélange le plus proche des cibles possibles à partir de l'optimiseur. Cet 
optimiseur reçoit la composition du mélange en entrée. L'optimiseur permet 
d'optimiser les proportions des courroies à partir d'un algorithme de contrôle. 


Cet algorithme de contrôle comprend une étape d'estimation des compositions 
des courroies et du mélange au prochain pas de temps et une étape de 
minimisation de l'écart entre le mélange estimé et les cibles an d'obtenir les 39 
proportions des courroies optimales 
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Les cibles pour 
la composition 
du mélange au cru sont présentées au tableau suivant . Elles sont 
représentatives des cibles utilisées dans certaines usines du groupe Lafarge. 























CaO 43,0 
SiO2 13,0 
A1203 3,9 
Fe203 1,3 





L'erreur d'estimation est calculée en faisant la moyenne des erreurs carrées 
entre la valeur estimée de la proportion d'oxyde dans la courroie ou le mélange 
et la vraie valeur de la proportion d'oxyde dans la courroie ou le mélange. 
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où i représente le numéro de la courroie, m le mélange ox, l'oxyde contrôlé et 
nox le nombre d'oxydes contrôlés 


Une stratégie pouvant être utilisée est de garder les proportions des courroies 
fixes. Cette stratégie peut être qualifiée d'absence de contrôle. Dans ce cas, les 
proportions des courroies qui restent fixes durant toute la simulation sont 
calculées en minimisant des équations des moyennes des compositions des 
courroies qui sont connues. Le tableau suivant présente les proportions de 


courroie. 





Argile 


0,0603 





Calcaire silicieux 


0,1889 





Calcaire 


0,7461 








Fer 








0,0048 





En absence de contrôle, le signal du mélange à la sortie est présenté par 
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(a) % d'oxyde de calcium dans le mélange 
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(c) % d'oxyde d'aluminium dans le mélange 


10- Fonction coût : 
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(d) % d'oxyde de fer dans le mélange 


Le modèle du procédé permet de prédire la composition du mélange à la sortie 
du broyeur et à la sortie du silo d'homogénéisation, le cas échéant. Ces valeurs 
prédites permettent le contrôle grâce à la minimisation d'une fonction de coût. 
Par exemple, chez Kural et Ozsoy (2004), la fonction de coût est dénie comme 


suit J = (ypr — yr) T (ypr — yr) 


où ypr est un vecteur représentant la composition du cru prédit et yr est le 
vecteur représentant les cibles Tandis que dans Ozsoy et al. (2002), la fonction 
de coût permet de prendre aussi en compte les coûts (ou la consommation 
d'énergie) des matériaux utilisés dans le mélange. Pour Lin et al. (1991) la 
fonction de coût à minimiser est très similaire, mais ce sont les coûts des 
matériaux utilisés qui sont en jeu plutôt que la consommation d'énergie. De 
plus, un tableau heuristique qui permet à l'opérateur de modifier la 
pondération dans la fonction de coût a été développé à partir de résultats de 
simulations qui ont montré que, lors de certaines perturbations, il est 
préférable de faire varier l'importance des différentes variables dans la fonction 
de coût. L'ajout de ce système auxiliaire a augmenté la robustesse du système 
(Lin et al., 1991). La fonction de coût inclut exclusivement les coûts reliés aux 
matériaux dans le système développé par Iino et al. (2008) Par contre, trois 
groupes de contraintes sont appliquées sur le mélange an que celui-ci atteigne 
les cibles qui sont exprimées en modules de fer, de silice et module 
hydraulique. Le premier groupe de contraintes lie les modules aux quatre 
oxydes qui composent le mélange, le deuxième groupe contraint les 
proportions des quatre courroies de sommer à 100% et le dernier groupe 
contient l'équation qui lie la composition du mélange aux proportions des 
courroies. 


11.Conclusion : 


Le système de contrôle présentement en place dans les cimenteries Lafarge 
(QMC) semble permettre un meilleur contrôle système de contrôle avec 
estimation de la composition du mélange et des matériaux par CKO ou CKS). 
Par contre, de petites modifications doivent y être apportés. Premièrement, il a 
été observé qu'effectuer le contrôle sur autant de cibles que de matériaux 
arrivant au mélange a un effet bénéfique par rapport à un contrôle effectué sur 
une cible de moins que le nombre de matériaux arrivant au mélange (ce qui 
est présentement le cas dans les cimenteries Lafarge). Deuxièmement, un 
paramètre de gain devra être rajouté sur le contrôle an d'amortir les réponses 
aux fluctuations. Enfin, les seuils sur les proportions des courroies devront être 
ajustées en fonction des proportions moyennes que les courroies risquent 
d'atteindre. 

Par ailleurs, aucun type de contrôle n'a semblé résoudre le problème des 
variations brusques d'un matériau arrivant au mélange. L'utilisation de 
moyennes mobiles n'a pas, non plus, permis de résoudre ce problème. Il est 
donc possible d'envisager certains travaux futurs. Premièrement, il serait 


intéressant d'effectuer des essais de contrôle sur un mélange de ciment réel. Il 
faudra déterminer si le choix d'un gain à partir de simulations des matériaux 
arrivant au cru est envisageable et si le contrôle avec un gain, ainsi choisi, 
améliore l'atteinte des cibles. 
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